Joint estimation of cosmological parameters from CMB and IRAS data 
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ABSTRACT 

Observations of large scale structure (LSS) and the Cosmic Microwave Background 
(CMB) each place separate constraints on the values of cosmological parameters. We 
calculate a joint likelihood based on various CMB experiments and the IRAS 1.2Jy 
galaxy redshift survey and use this to find an overall optimum with respect to the free 
parameters. Our formulation self-consistently takes account of the underlying mass 
distribution, which affects both the CMB potential fluctuations and the IRAS redshift 
distortion. This not only allows more accurate parameter estimation, but also removes 
the parameter degeneracy which handicaps calculations based on either approach alone. 
The family of Cold Dark Matter (CDM) models analysed corresponds to a spatially- flat 
universe with an initially scale- invariant spectrum and a cosmological constant. Free 
parameters in the joint model are the mass density due to all matter (fJ m ), Hubble's 
parameter (h = Hq/IOQ km s _1 Mpc~ 1 ), the quadrupole normalisation of the CMB 
power spectrum (Q) in /iK, and the IRAS light-to-mass bias (bi ias ). Throughout 
the analysis, the baryonic density (f2&) is required is to satisfy the nucleosynthesis 
constraint tl^h 2 = 0.024. Results from the two data sets show good agreement, and the 
joint optimum lies at f2 m = 0.39, h = 0.53, Q = 16.96 fiK, and &i ras = 1.21. The 68 per 
cent confidence intervals are: 0.29 < Q m < 0.53, 0.39 < h < 0.58, 15.34 < Q < 17.60, 
and 0.98 < 

^iras < 1.56. For the above parameters the normalisation and shape of 
the mass power-spectrum are erg = 0.67 and T = 0.15, and the age of the Universe is 
16.5 Gyr. 
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1. Introduction 



Astronomical observations allow us to evaluate cosmological models, and determine likely 
values for the parameters within them. As the variety and depth of these observations have 
improved, so too have the techniques for comparing them against theoretical predictions. 

Observations of anisotropies in the Cosmic Microwave Background (CMB) provide one of the 
key constraints for cosmological models and a significant quantity of experimental data already 
exists. By comparing the power spectrum of the CMB fluctuations derived from these experiments 
with the power spectra predicted by different cosmological models it is possible to set constraints 
on the value of certain cosmological parameters (e.g. Hancock et al. 1998] , [Lineweaver et al. 1997 , 
and references therein). 

Galaxy redshift surveys, mapping large scale structure (LSS), provide another cosmologically 
important set of observations. The clustering of galaxies in redshift-space is systematically 
different from that in real-space ( [Kaiser 1987 ). The mapping between the two is a function of the 
underlying mass distribution, in which the galaxies are not only tracers, but also velocity test 
particles ( Lahav 1996 ). Many techniques have been developed for estimating this mapping ( |Yahil 
let al. 1991| , || aiser et al. 199l|) . Statistical quantities can be generated for a given cosmology and 
these used to constrain model parameters through comparison with survey data ( [Fisher, Scharf fc| 
Lahav 199j |Cole et al. 1995| , [Heavens k Taylor 1995| , [Fisher fc Nusser 1996| , |Willick et al. 1997| ). 



Estimates derived separately from each of these two data sets have problems with parameter 
degeneracy. In the analysis of LSS data, there is uncertainty as to how well the observed light 
distribution traces the underlying mass distribution. The light-to-mass linear bias, b, introduced to 
account for this uncertainty, affects the value of many central cosmological parameters, and makes 



any identified optimum degenerate (Strauss & Willick 1995). Similarly, on the basis of CMB data 
alone, there is considerable degeneracy ( [Bond et al. 1995 (j ) between h = H /1Q0 km s^Mpc" 1 
and the energy density Q\ due to the cosmological constant ( |Carroll, Press &: Turner 1992 ). This 
leads to poor estimation of the baryon (fib) and total mass (f2 m ) densities ( [Lineweaver et al. 1997 ). 



In this letter, we combine results from a bandpower approach covering a range of CMB 
experiments, with a likelihood analysis of the IRAS 1.2Jy survey, performed in spherical harmonics. 
We present a self-consistent formulation of CMB and LSS parameter estimation. In particular, 
our method expresses the effects of the underlying mass distribution on both the CMB potential 
fluctuations and the IRAS redshift distortion. This breaks the degeneracy inherent in an isolated 
analysis of either data set, and places tight constraints on several cosmological parameters. Indeed, 
it is unsurprising that the two data sets are complementary, given that they sample our universe at 
extreme ends of its evolution. For simplicity, we restrict our attention to inflationary, Cold Dark 
Matter (CDM) models, assuming a flat universe with linear, scale-independent biasing. Other 
recent studies which combine CMB and LSS include Gawiser & Silk 1998 and [Eisenstein, Hu fc 



Tegmark 1998 . 



- 3 - 



2. CMB Parameter Estimation 



2.1. Experimental Data 



Since the discovery of CMB fluctuations by the COBE satellite ( [Smoot et al. 1992 , Bennett 
et al. 1996| ), several other experiments have also measured CMB anisotropics over a wide range 
of angular scales. These experiments include ground-based beam-switching experiments such 
as Tenerife ( |Hancock et al. 1997| , [Hancock et al. 1994| ), Python ( [Ruhl et al. 1995| ), South Pole 
( Gundersen et al. 1995| ) and Saskatoon (Netterfield et al. 1997); balloon-borne instruments 
such as ARGO ( pe Bernardis et al. 1994[ ), MAX ( |Tanaka et al. 1996| ), and MS AM QCheng et 
al. 1994, Cheng et al. 199"6| ) and the ground-based interferometers CAT ( [Scott et al. 1996| ) and 
OVRO QLeitch et al.l99S| ). 



These observations have resulted in a first estimate of the CMB power spectrum and they are 



discussed in more detail by [Hancock et al. 1998| . In particular, Hancock et al. display the window 
function W\ for each experiment and convert the level of anisotropy observed in each case to flat 
bandpower estimates (A7]/T) ± a centered on the effective multipole Z e ff of the corresponding 
window function (see below). The resulting CMB data points are plotted in Fig. |l|, together with 
their 68 per cent confidence limits. These confidence limits have been obtained using likelihood 



analyses and hence incorporate uncertainties due to random errors, sampling variance (Scott 



ig variance ( 


Scott, 


Scaramela & 


Vittorio 



1993). A discussion of possible additional uncertainties due to contamination by foreground 
emission is given by Rocha et al. 199S| , 

The data points plotted in Fig. [I] differ slightly from those given in Hancock et al. 1998 as 
follows. The old Python point has been replaced by two new points corresponding to the Py Ills' 
and Py I, Py II and Py IIIl observations respectively ( Piatt et al. 1997| ). The Tenerife point has 
been updated following the analysis of the full two-dimensional data-set and a detailed treatment 
of atmospheric effects ( Gutierrez 1997| ). The Saskatoon points have been increased by 5 per 
cent, as suggested by the recent investigation of systematic calibration errors in the Saskatoon 
experiment (Leitch, private communication). We have added the results from the second CAT 
field, reported in [Baker et al. 1998] and also the OVRO point reported by Leitch et al. 1998. 



2.2. Method 



Temperature fluctuations in the CMB are usually described in terms of the spherical harmonic 
expansion 



oo I 
1=0 m=-l 



(1) 



- 4 - 



from which we define the ensemble-average CMB power spectrum C\ = (|a; m | 2 ). Alternatively, we 
may describe the fluctuations in terms of their autocorrelation function 

The power in the CMB fluctuations observed by an experiment with window function W; is then 
given by 

c obs (o) = ^ (^) 2 ) = ^ £ (« + 1) c i wi , (3) 

and for each experiment we define the flat bandpower by 



T V /(Wi) ' 
where J(Wj) is defined (|Bond 1995 oj |Bond 1995 1 as 



AT,_ /C obs (0) 



= y iil^ ^ . (5) 

This flat bandpower estimate is centred on the effective multipole 

We wish to compare these bandpower estimates with those predicted by different cosmological 
models. Varying the values of model parameters, we calculate corresponding, predicted CMB 



power spectrum Cj, using the Boltzmann code of Seljak Zaldarriaga 1996 , This is then used to 
calculate the predicted flat bandpower A/] for each experiment. The chi-squared statistic for a 
given set of parameter values, <3 cm b, is then 



X 2 («cmb) = E ^([ AT ° bS l " [ATP red (a cmb )j J , (7) 



where N<i is the number of CMB data points plotted in Fig. |l| (20 in this analysis). Moreover, 
since the CMB data points plotted in Fig. |l] were chosen such that no two bandpower estimates 
come from experiments which observed overlapping patches of sky and had overlapping window 
functions, we may consider them as independent estimates of the CMB power spectrum. As the 
cosmic variance has already been taken into account in deriving the flat bandpower estimates, the 
likelihood function is given simply by £ cm b oc e~^ I 2 . 

As mentioned above, we assume that the Universe is spatially flat, and that there are no 
tensor contributions to the CMB power spectrum. We take the primordial scalar perturbations 
to be described by the Harrison-Zel'dovich power spectrum for which n s = 1, and further assume 
that the optical depth to the last scattering surface is zero. 
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The normalisation of the CMB power spectrum is determined by Q, which gives the strength 
of the quadrupole in /j,K, such that 

4vr / Q x 2 



where To is the average CMB temperature. The expansion rate of the Universe is given by 
Hubble's parameter h = km s -1 Mpc , while fi cc i m an d fib denote the density of the 

Universe in CDM and the baryons respectively, each in units of the critical density. Given that we 
assume a flat universe, but investigate models where 

fim = (ficdm + fi b ) < 1 , (9) 
the shortfall is made up through a non-zero cosmological constant A such that 



fi A = i - fi m = — . (io) 

orir. 



A 

Furthermore, we restrict our attention to models that satisfy the nucleosynthesis constraint 
fi^/i 2 = 0.024 ( Tytler, Fan fc Buries 1996 ). Thus we consider the reduced set of CMB parameters 

«cmb = {Q,h,VL cdm } . (11) 

In Section 4, we derive a set of joint parameters linking these to the IRAS parameter set. 



3. IRAS Parameter Estimation 



3.1. IRAS 1.2Jy Survey 



The IRAS surveys are uniform and complete down to Galactic latitudes as low as ±5° from 
the Galactic plane. This makes them ideal for estimating whole-sky density and velocity fields. 
Here, we use the 1.2 Jy IRAS survey (Fisher et al. 1995), consisting of 5313 galaxies, covering 
87.6% of the sky with the incomplete regions being dominated by the 8.7% of the sky with \b\ < 5°. 

In principle, the method we are using can be extended to account explicitly for the incomplete 
sky coverage. However, we adopt the simpler approach of smoothly interpolating the redshift 
distribution over the missing areas ( Yahil et al. 1991~1) . The effects of this interpolation on the 



computed harmonics have been shown to be negligible (Lahav et al. 1994). 



3.2. Method 

In this letter, we assume linear, scale-independent biasing, where &i ras measures the ratio 
between fluctuations in the IRAS galaxy distribution and the underlying mass density field: 

(<Wp)iras = & iras (<Wp)m ■ ( 12 ) 
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We note that biasing may be non-linear, stochastic, non-local, scale-dependent, epoch-dependent 
and type-dependent (e.g. Dekel and Lahav 1998| , |Tegmark fc Peebles 1998 , Pen 1998] , [Bagla 
1998, Blanton et al. 1998, and references therein). For a linear bias parameter, 



the 



velocity and density fields in linear theory (Peebles 1980) are linked by a proportionality factor 
/3 iras = r2 m °' 6 /friras, such that 

V-V = -F Aras (S P /p) m (13) 

Statistically, the fluctuations in the real-space galaxy distribution can be described by a 
power spectrum, P^(k), which is determined by the rms variance in the observed galaxy field, 
measured for an 8 /i _1 Mpc radius sphere (o"8,iras) and a shape parameter (e.g. T in equation IS). 
The observed 08,iras is related to the underlying as for mass through the bias parameter, such that 

C8,iras = &iras 0"8- 

The approach we use in this letter follows Fisher, Scharf & Lahav (1994; hereafter FSL), and 
we include here only a brief introduction to our technique. FSL provides a detailed description of 
the spherical harmonic approach to parameter estimation. A flux-limited, redshift-space density 
field can be decomposed into spherical harmonics with coefficients 

v„ 



a ir, 



J2f(si)Y lm (§j) 



(14) 



i=i 



where 2V g is the number of galaxies in the survey, and f(s) is an arbitrary radial weighting 
function — this process is analogous to Fourier decomposition, but instead using spherical basis 
functions. The sum over galaxies in this equation can be rewritten as a continuous integral of the 
density fluctuation field S(s) over redshift-space: 



Hm 



d 3 s <f>{r)f{s) [l + £ s (s)] Y lm (s) 



(15) 



where 4>(r) is the radial selection function of the survey, evaluated at the real-space distance of the 
i th galaxy. 

As detailed in FSL, assuming the perturbations introduced by peculiar velocities are small, 
the expected linear theory values for the harmonic coefficients are 

oo 



a l, 



= - I dkk 2 P K {k) *f (fc) + Arastfftfc) 

LT w 



(16) 



Here, X I ,R is the real-space window function, while Aras^ is a "correction" term which embodies 
the redshift distortion. 

For a given set of cosmological parameters, redshift-space harmonics aj can be calculated 
for different weighting functions f J (r), j = 1, ...,N. If the underlying density field is Gaussian, on 



the basis of the coefficients predicted in equation 16, the likelihood of the survey harmonics can 
be calculated as 

£ iras oc |A|"2 exp (-- [a T A" 1 a] J . (17) 
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Here a is the vector of observed harmonics for different shells and A is the corresponding 
covariance matrix, which depends on the predicted harmonics. In addition to the (|of m |)LT f rom 



equation 16, these predicted harmonics will also have a shot-noise contribution (|a; m |) SN due to 



the discreteness of survey galaxies. Note that the argument of the exponent in equation 17 is 
simply (— x 2 /2), and that here the normalisation of the likelihood function does depend on the 
free parameters (unlike in the CMB likelihood function). 

In this letter, we follow FSL and use four Gaussian windows centered at 38, 58, 78, and 98 
/i _1 Mpc each with a dispersion of 8 ft. Mpc. For each window, we compute the corresponding 
weighted redshift harmonics for the IRAS 1.2Jy catalog, and use these to determine the likelihood 
of a given set of parameter values. Since our analysis is valid only in the linear regime, we restrict 
the likelihood computation to Z max = 10 (corresponding to 120 degrees of freedom). Hence, the 
IRAS likelihood function has a parameter vector 

"iras = {Aras, 0"8,iras, T} . (18) 

Again, the linkage between these and the CMB parameters is discussed in section 4 below. 



4. Joint analysis 



Given the large number of parameters available between the two models, it is important both 
to find links for joint optimisation, and to decide which parameters can be frozen. From section 3, 
we have six variables between the two models: {Q, h, O cc j m , Aras> C8,iras> T}. These can be reduced 
further by expression in terms of core cosmological parameters. The IRAS normalisation can be 
calculated as as = f{Q m ,Q,T) flBardeen et al. 1986; , jEfetathiou, Bond fc White 1995Q , while the 
CDM shape parameter is well approximated flSugiyama 1991) by 



r = O m h exp —Ob 



1 + 



O r , 



Meanwhile, we have shown above that O m 
Hence, the final, joint parameter space is 



O c dm + Ob, Aras — &m /&irasj an d 0"8,i) 



C^joint — {hi Qi Orru ^iras} • 



(19) 

08&iras- 

(20) 



As the IRAS and CMB probe very different scales and hence are assumed to be uncorrelated, the 
joint likelihood is given by 

In (£joint) = In (£cmb) + In (Aras) • (21) 



5. Results 

The complementary nature of the two data sets is demonstrated in Fig. which shows 
likelihood contours in the {O m , /i}-plane after marginalising over Q and b\ ras . The fundamental 
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CMB-side degeneracy in {Cl m , h} is seen in the flat trough running across Fig. |2|a. The IRAS 
degeneracy is in a different direction (Fig. Qb) and the two data sets agree well in the region where 
the lines of degeneracy overlap. Combining the two data sets breaks the degeneracy and leads to 
a well defined joint optimum (Fig. |2]c) . The joint CMB plus IRAS likelihood function in the {fl m , 
Q}-plane and {h, Q}-plane is plotted in Figs || & |] respectively. In each case, a marginalisation 
has been performed over the other two parameters. 

The joint likelihood (equation |2l| ) was maximized with respect to the 4 free parameters 



(equation 20) using standard optimisation techniques ( Press et al. 1992 ) and the best fit parameters 
are shown in Table |]. For this set of parameters, we find the values of the reduced x 2 for the IRAS 
and CMB data respectively to be 1.18 and 1.03, confirming that both data-sets agree well with the 
models used. Taking the CMB and IRAS data together the total reduced x 2 is 1.16. Recalculating 



the joint optimum using the simpler formula T = O m h in place of that in equation 19 had little 



effect. Further, the optimum was robust to changes in IRAS / max in the range 4 < Z max < 10. 

To obtain 68 per cent confidence limits on each of the free parameters it is necessary to 
marginalise over the remaining free parameters. To achieve this we evaluated the joint likelihood 
function on a 4-dimensional grid of parameter values. The range of values and number of grid 
points in each direction were 5 < Q < 30, 50 steps; 0.3 < h < 0.9, 50 steps; 0.1 < $7 m < 1.0, 
100 steps; 0.7 < 6j ras < 2.0, 20 steps. For each parameter, the corresponding one-dimensional 
marginalised probability distribution was calculated by integrating the likelihood function over the 
other variables. The marginalised distribution for each parameter is shown in Fig. || in which the 
dashed vertical lines denote the 68% confidence limits quoted in Table ||. In general, the peak of 
the one-dimensional probability distribution for each variable will differ from the global optimum 
across all parameters. However, for all four variables in this system, the two values are found to 
be extremely close. 

In addition we evaluated the covariance matrix at the joint optimum. The covariance 
matrix is simply the inverse of the Hessian at the joint optimum. The Hessian is given by 
<9 2 ln (Cjoint) / (daidaj) for pairs of parameters Oj and ay and is evaluated using a standard 
central-difference algorithm. Taking the square-root of the diagonal elements of the covariance 
matrix, the standard errors on the parameters are found to be AQ = 0.67 fiK, Ah = 0.03, 
Af2 m = 0.04 and A6j ras = 0.1. Comparing these with the marginalised errors quoted in Table [l] we 
note that the marginalised errors are consistently larger, as expected. By rescaling the covariance 
matrix so that its diagonal elements equal unity, we obtain the correlation matrix shown in 
Table 0. We see that the most strongly correlated parameters are fi m and h. 



6. Discussion 



The results of this joint optimisation are in reasonable agreement with other current estimates. 
The relatively low value of f2 m ~ 0.4 is close to that found by others ( |White et al. 1993 , Bahcall, 
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0.29 < n m < 0.53 
0.39 < h < 0.58 
15.34 < Q < 17.60 
0.98<6 iras < 1.56 

Derived parameters 



n b 0.085 

o- 8 0.67 

0"8,iras 0.81 

r 0.15 

Aras 0.47 

Age (Gyr) 16.5 



Free parameters 
Q m 0.39 
h 0.53 
Q (/xK) 16.95 

^iras 1.21 



Table 1: Parameter values at the joint optimum. For the free parameters the 68% confidence limits 
are shown, calculated for each parameter by marginalising the likelihood over the other variables. 





Q 


h 




&iras 


Q 


1.00 


0.40 


-0.46 


-0.49 


h 


0.40 


1.00 


-0.82 


-0.57 




-0.46 


-0.82 


1.00 


0.09 


&iras 


-0.49 


-0.57 


0.09 


1.00 



Table 2: Parameter correlation matrix at the joint optimum 
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Fan and Cen 1997), and is in line with recent supernovae results ( Perlmutter et al. 1998j ). However, 
given the assumption of a flat universe, it requires a very high cosmological constant (p,\ = 0.6). 
Gravitational lensing measurements have constrained < 0.7 ( Kochanek 1996] , Falco, Kochanek] 
fc Munoz 1998| ). Our value for the Hubble constant, h = 0.53, agrees well with several other 
measurements flSugiyama 1995 , Lineweaver et al. 1997), but falls at the low end of the generally 
accepted range from local measurements ( [Freedman et al. 1994 ). Assuming the nucleosynthesis 
constraint n b h 2 = 0.024 flTytler, Fan k Buries 1996| , [Steigman, Hata & Felten 1998[) , the optimal 
baryon density is found to be S\ = 0.085. Our value for the combination <78^ m °' 6 = 0.38 is 
lower than the one derived from measurements from the peculiar velocity field, cxs^m ' 6 ~ 0.8 
( Frcudling et al. 1998| ). Our values are closer to the combination derived from cluster abundance 
<7 8 ft m ' 5 ~ 0.5 ( [Eke, Cole fc Frenk 1996| , |Bahcall, Fan and Cen 19971) . Finally, for spatially-flat 
universes the time since the Big Bang for the values of our Q m and h at the joint optimum is 16.5 
Gyr. 



On the IRAS side, Aras = 0.47 is in agreement with several other measurements ( [Strauss 



1989, [Schlcgcl 1995 , Willick et al. 1997), although there are other measurements which place Ai 



much higher (Dekel et al. 1993, Sigad et al. 1998). [Willick et al. 1997 discuss the discrepancies 
between the various measurement techniques, and why they lead to such distinct results. Finally, 
the IRAS mass-to-light bias is seen to be slightly greater than unity (&i ra s = 1-21), suggesting 
that the IRAS galaxies (mainly spirals) are reasonable (but not perfect) tracers of the underlying 
mass distribution. We also note that our joint IRAS+CMB optimal values for o^irasi Aras and 
r (Table 1) are in perfect agreement with the values derived from IRAS alone (FSL, [Fisher 
|1994[ ). However at fixed 08,iras = 0.69 based on the IRAS correlation function FSL found a higher 
Aras = 0.94 ± 0.17 and T = 0.17 ± 0.05(1 — a). We emphasise that the naive linear biasing should 
be generalised to more realistic scenarios (e.g. pekel and Lahav 1998 ). 

As discussed in Section |j| the total CMB+IRAS reduced x 2 at the joint optimum is 1.16, 
indicating that the optimum model is a good fit to both data-sets. We may compare this value 
with that obtained by Gawiser fc Silk 1998j . Using several data-sets Gawiser &: Silk 1998 found the 
reduced x 2 for ACDM to be 1.9, as opposed to a value of 1.2 for their 'best' model of Cold+Hot 
Dark matter (CHDM). However, using only the CMB and IRAS x 2 values for ACDM in their 
Table 3, the total reduced \ 2 ls found to be 1.00 as compared to a value of 0.95 for CHDM. Thus, 
using only the CMB and IRAS data-sets, ACDM is found to fit the observations as well as CHDM, 
and the results quoted by pawiser fc Silk 1998 are consistent with those presented here. 

The near future will see a dramatic increase in LSS data (e.g. the PSCZ, SDSS, 2dF surveys) 
and detailed measurements of the CMB fluctuations on sub-degree scales (e.g. from the Planck 
Surveyor and MAP satellites). These will allow more accurate parameter estimation and the 
untying of the various parameters held fixed in the present work. 
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Fig. 1. — The data points used in the calculation of the CMB likelihood function. The overlaid 
curve is a model evaluated with optimal parameters shown in Table 1. 
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Fig. 2. — The likelihood function in the {fi m , /i}-plane, after marginalisation over Q and &i ras , for 
(a) CMB alone, (b) IRAS alone and (c) joint CMB and IRAS. The contours denote the 68, 90, 95 
and 99 per cent confidence regions. 
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Fig. 3. — The joint CMB+IRAS likelihood function in the {fi m , Q}-plane after marginalisation 
over h and &i ras . The contours denote the 68, 90, 95 and 99 per cent confidence regions. 
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Fig. 4. — The joint CMB+IRAS likelihood function in the {h, Q}-plane after marginalisation over 
U m and 6i ras . The contours denote the 68, 90, 95 and 99 per cent confidence regions. 
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Fig. 5. — The one-dimensional marginalised probability distributions for each of the four 
parameters. The vertical dashed lines denote the 68% confidence limits. The horizontal plot 
limits are at the 99% confidence limits. 



